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Abstract. Quantum models on the hyper-cubic d-dimensional lattice of spin-i particles 
interacting with linear oscillators are shown to have three ferromagnetic ground state order 
parameters. Two order parameters coincide with the magnetization in the first and third 
directions and the third one is a magnetization in a continuous oscillator variable. The 
proofs use a generalized Pcicrls argument and two Griffiths inequalities. The class of spin- 
oscillator Hamiltonians considered manifest maximal ordering in their ground states. The 
models have relevance for hydrogen-bond ferroelectrics. The simplest of these is proven to 
have a unique Gibbsian ground state. 
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1 Introduction 

In this paper we consider quantum lattice models of oscillators interacting with spins whose vari- 
ables are indexed by the sites of a hyper-cube A with the finite number of sites |A| in the 
hyper-cubic lattice 7L d . Interaction is considered to be short-range and translation invariant. 
The corresponding Hamiltonian H\ is expressed in terms of the oscillators variables q\ = 
(q x ,x 6 A) € K' A ' and spin \ Pauli matrices S l A = (S l x ,x G A, I = 1,3), defined in the ten- 
sor product of the 2l A l-dimensional Euclidean space and the space of square integrable functions 
h\ = (®E 2 )I A I ®L 2 (IRl A l), 

H A = J2[-d 2 x +^ 2 (q x + riMSl)) 2 -^} + ^2jASl A] +V A , ft > 0, r, € M, (1.1) 

xeA aca 

where d x is the partial derivative in q x , J a and V\ are real- valued measurable functions, the 
first of which depends on q\ and the second on S\, q\ and the translation invariant 4> x is given 
by (see Remark [2] in the end of the paper) 

<Px{sa) = ^2 Mx;A)s [A] , J (x;x) = l. (1.2) 
ACA 

*This paper is a contribution to the Proceedings of the Seventh International Conference "Symmetry in 
Nonlinear Mathematical Physics" (June 24-30, 2007, Kyiv, Ukraine). The full collection is available at 
|http://www.emis.de/journals/SIGMA/symmetry2007.html| 
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For products of operators, functions and variables we use the following notation: Bui = Y\ B x . 

The scalar product in ((g>E 2 )l A l, L, A will be denoted by (-,-)oj (v)> respectively. The Schwartz 
space of test functions on M. n will be denoted by S(IR n ). 

We require that the Hamiltonian is well defined and bounded from below on (£g)E 2 )l A l (g) 
Co°(]R' a I), i.e. the tensor product of the 2' A '-dimensional Euclidean space and the space of 
infinitely differentiable functions with compact supports. The ground state average for an 
observable (operator) F is determined by 

(F) A = zxH*a,f*a), z a = (\P a ,*a) = II^aII 2 , 

where Z\ is a partition function. For partial cases of F we have 

(Q[A]<P[A'](Sa))a = Z^ 1 J q[A](^A(qA),^[A'](Sl)^A(qA))odqA, 

where the integration is performed over M' A ' and q\A\ is the operator of multiplication by q x . 

We will employ the orthonormal basis ^9(sa) of the Euclidean space (<8>1E 2 )' A ', diagonali- 
zing S A , which is chosen in the following way: i/j a (s\) = l S' X £A' l l j0 (s x ), s x = ±1> V^U) = (1)0), 
</>°(-l) = (0,1), 5V°(s) = *l>°(s), S 3 4>°{s) = sip°(s). For F € L A we have the following 
decomposition 

%) = J%;«aK(«a), 

SA 

where the summation is performed over the |A|-fold Cartesian product (—1, l)' A 'of the set (—1, 1). 

The Hamiltonians in (jl.ip are employed in hydrogen-bond ferroelectric crystal models, consi- 
dered in [H El E] , and describe interaction between heavy ions (oscillators with constant fre- 
quency) and protons (spins). The second term with J a = 0, \A\ > 2 corresponds to the energy 
of protons, tunneling along hydrogen bonds from one well to another, and J x is associated with 

the tunneling frequency. The last term in the expression for our Hamiltonian Va = £ Ji(A)Sf A] 

ACA 1 1 

describes many-body interaction between protons [J\(A) is the intensity of the |j4|-body inter- 
action) . 

A rigorous analysis of a mean-field version of the Hamiltonian in (jl.ip with <p x {S\) given 
in CL2]) that is linear in S 3 , J x ^ 0, V A = -(w) 2 £ 4>l(S A ), J a = for \A\ > 2 and J {x;y) 

x£A 

uniformly in lattice sites proportional to |A| _1 was carried out in [3] in the framework of the Bo- 

golyubov approximating Hamiltonian method [J] and occurrence of spin and oscillator orderings 

(the corresponding order parameters are non-zero) for two-body interaction between protons at 

non-zero temperatures was proved. To establish such orderings for ground states without the 

mean-field limit in a general case is an important task for a theory. 

The oscillator and spin orderings are established if one proves the existence of ferromagnetic 

oscillator and spin long-range orders (lro's) in ground states for the corresponding Hamiltonians. 

This means that the ground state averages {q x q y )A, {S x Sy)\, j = 1, 3 are bounded uniformly in A 

from below by positive numbers . Occurrence of the ferromagnetic lro's implies the existence of 

the spin order parameters (magnetizations in the first and third directions) M l K = |A| _1 £ S l x , 

xeA 

I = 1,3, and the oscillator order parameter Q\ = | | 1 £ q x in the thermodynamic limit 

xeA 

(A — ► Z d ) since the ground state averages of their squares are uniformly bounded in A from 
below by a positive number. 
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In this paper we find functions V A , depending on oscillator variables, for which ground states 
or eigenstates ^a of the Hamiltonians in are Gibbsian 

*a(<?a) = £V ^ 8a « a ^(*a)^M?a), (1.3) 

with the linear in ^a spin-oscillator quasi-potential energy U 

U(s A ; q A ) = 2(iri q x <t) x {s A ) + U°(s A ), (1.4) 

where iPoa(qa) = I ipo(q x ), ipo(q) = (m 71 " -1 ) 3 exp{ — ^q 2 }, is the ground state of the free 

oscillator Hamiltonian, that is the first term in the right-hand side of with r/ = 0. We 

prove the maximal ordering in the corresponding systems (provided some simple conditions 
on U° and 4> x are satisfied): magnetizations Mi, I = 1,3, Q A are non-zero. No other ground 
states are known with such the property. 

Gibbsian ground states were introduced by Kirkwood and Thomas in [5] in XZ spin-^ models 
with Hamiltonians (linear in 5 1 ) that include the spin part of (jl.ip . i.e. the second and third 
terms, with only J x ^ and the periodic boundary condition (this boundary condition is not 
essential). Matsui in O E] enlarged a class of spin-1 XZ-type models in which Gibbsian states 
exist. The method was further developed by Datta and Kennedy in [Sj. An application of 
classical spins systems for constructing of quantum states was given in [9]. 

In [10] we showed how to find V A for a given Gibbsian ground state and established existence 
of lro's in S 1 and S 3 for a wide class of the spin-^ XZ models (see also [HUTS]). This reference 
contains the most simple proofs of the existence of lro's in ground states of quantum many-body 
systems. A reader may find a review of the results concerning several quantum orders in it (see 

also pa Hang). 

The ground state in (jl .31) can be represented in the following equivalent form 

^A(gA) = ^e-^*^)<( SA )V'0A(gA + #A(sA)), U, = U° - /V £ (1.5) 

sa xeA 

From (jl.5|) . orthonormality of the basis (see the beginning of the third section) and the 
equalities J ipQ(q)dq = 1,J ipQ(q)qdq = it follows that 

{qxq y )A = ^ 2 (0x(5'a)^?/('S'a))a = V 2 (4>x{vA)4>y(vA))*A, (S x Sy) A = {a x a y )^ (1.6) 
where cr x (s A ) = s x and 

(M°a)M°a))*a = z;lY, e ~ u * (SA) ^ s ^M s A), z* A = ^e- u * (SA) . 

SA S A 

Equalities in (jl.6p reduce a calculation of averages in our quantum systems to a calculation of 
averages indexed by a star in Ising models. For a short-range ferromagnetic potential energy U* 
the first Griffiths inequality holds: (<J[a])*a > 0. As a result the following statement (principle) 
is true. 

Proposition 1. Let Jq > in (|1.2p and ferromagnetic Iro occur in the Ising model with the 
ferromagnetic potential energy U* given by (jl.5p . that is (a x <Jy)*\ > uniformly in A, then 
ferromagnetic Iro occurs in oscillator variables and S 3 in the quantum spin-oscillator system 
with the Hamiltonian (jl.ip and ground state ^/ A in (j 1 . 3 j) . 

Usefulness of Gibbsian ground states is explained by comparative simplicity of a proof of 
existence of Iro. Our results show that Gibbsian ground states are expected to appear in many 
quantum spin-oscillator systems with non-trivial interactions. 

Our paper is organized as follows. In Section [5] we formulate our main results in a lemma 
and theorem. In next sections we prove them. 
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2 Main result 

We establish that Gibbsian ground states exist for J a depending on q\ if 

Va = ~ JAe-3 WA &\ W A (Sl) = U(Sl A ;q A ) - U(Sl;q A ), (2.1) 
ACA 

where S\ A = (—S A ,S A \ A ). V A determines an unbounded operator. Negative J a generate 
positive functions V A and this enables us to prove the following lemma. 

Lemma 1. Let V A be given by ()2.ip and have a domain D(V A ), J a be bounded negative functions 
and U in (|2.ip coincide with U in (|1.4p Then H A is positive definite and essentially self-adjoint 
on the set ((g) C 2 ) |A| <g>S(Rl A l) HD(V A ) that contains (<g)C 2 ) |A| <g>C£°(M |A| ) and^ A , given by (TOll . 
is an eigenfunction of the Hamiltonian (jl.ip with the zero eigenvalue and is its (unique) ground 
state if J a < (if the uniform bound J x < J_ < holds). 

Remark 1. If the functions J a are only negative then ^f A is the ground state of the self- adjoint 
extension of the Hamiltonian preserving positive definiteness. 

The most simple translation invariant short-range U° is ferromagnetic 

C/° (sA) = _ a ^jO S[A]) j0> 0i (2 . 2) 

ACA 

where J A = for odd \A\, a > 0. 

Theorem 1. Let all the conditions of Lemma [1] be satisfied, Jq > 0, Jo(x;A) = for even \A\, 
J (0; 1), Jq ! > J and U° be given by ([221). Then 

I. For a sufficiently large (3 = {rfn + a) J > 1 there exist ground state Iro's in and oscillator 
variables for d>2; 

II. Let the positive constants C, Bj, j = 1,2, independent of A, exist such that \<fi x (s A )\ < C 

w ( /\ SA ) = \€>(4) - A(«a)| < B h j = 1, 2, 

x'eA 

where \A\ = 2. Then ground state ferromagnetic Iro in S 1 occurs in arbitrary dimension d. 
If U° = then V A from (|2.ip is given by 

Va = - ^2 J AV[A], v x = coshu x + S 3 smliu^, u x = 2rjn4> x (q A ). 

ACA 

This equality follows from the equalities 

W A (Sl) = -2^2u x S'i e aS3 = cosha + 5 3 sinha, (S 3 ) 2 = /. 

In the simplest case the conditions of item II of Theorem Q] can be checked without difficulty. 
Proposition 2. Let cp x be linear in <S 3 in ()1.2|) and ||Jo||i = ^2 1^0(^)1 < oo, where the sum- 

X 

motion is performed over 7L d . Then the conditions of item II of Theorem^ are satisfied. 
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Note that if one uses the Pauli matrices with instead of the unity as matrix elements 
then V A should be changed by adding to Wa the number — \A\ In 2. 
If one chooses the anti- ferromagnetic U° in (|2.ip . specifically, 

U°(s A ) = fir] 2 0x( s a) + a s x s y , a > 0, 

xeA (x,y)eA 

where x, y are nearest neighbors, then it can be easily proved that the spin lro in the third 
direction will be anti-ferromagnetic, generating a staggered magnetization (spins at the even 
and odd sublattices take different values). 

The interesting and important property of the Hamiltonians with V A given by (|2.ip is that 
they are simply related to generators of stationary Markovian processes (see, also, |16|). We 
believe that it is possible to apply the same mathematical technique for proving existence of 
order and phase transitions in equilibrium quantum systems and non-equilibrium stochastic 
systems (see [17]), 

3 Proof of Lemma [I] 

For our purpose it is convenient to pass to a new representation. It is determined by the Hilbert 
space of sequences of functions F(q A ;s A ), s x = ±1, which are found in the expansion of the 
vector F £ L A mentioned at the beginning of the introduction, with the scalar product 

{Fi,F 2 ) = y~) / F 1 (q A ;s A )F 2 {q A ;s A )dq A , 

SA J 

(F 1 (q A ),F 2 (q A )) = ^Fi(g A ; s A )F 2 (q A ; s A ), 

SA 

where the integration is performed over Rl A l. Here we took into account the orthonormality of 
the basis, i.e. the equality 

«(sa), K(s' a )) = 5(s A ; s' A ) = J] S SxA , 

xeA 

where 5 StS > is the Kronecker delta. Let 

h x = -d 2 x + fi 2 (q x + V<Px(S A )) 2 ~ 

h °x = ~ d l + V 2 ql - V = (-9 X + mx){d x + fiq x ) 

and 

h A = y £h x , h° A = ^2h° x . 

xeA xsA 

(T x (4>)F) (q A ) = F (q x + r]^ x (s A ), q A \ x ;s A )ip A (s A ) , 

S A 

then 

h x = T x {(t>)h x T- 1 (<f>), h A = T A h\T^\ 

Here we took into account that differentiation commutes with T A . 
Our Hamiltonian is rewritten as follows 

H A = h A +Y^ JaP [A] , Pa = S{ A] - e-^ s l\ 

ACA 
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The remarkable fact is that the symmetric operator Pa and the harmonic operator h x have both 
common eigenvector ^a with the zero eigenvalue (see Remark [3] in the end of the paper) . Note 
that the space of ground states of the operator h x (eigenfunctions with the zero eigenvalue) is 
2l A l-fold degenerate since S% is diagonal and the Laplacian is translation invariant. From (|1.4p 
and the definition of T x ((j)) if follows that T~ l {(f))^ A is equal to 4>o{q x ) multiplied by a function 
independent of q x 

T-\4>)* A = XX^^A^A^OAteAV* +V<t>A\x{sA\x))M<lx)- 

Hence h^T' 1 (fyV A = and h x V A = T^h^T- 1 ^)^ A = 0. The proof that V A is an eigenvec- 
tor with the zero eigenvalue of Pa is inspired by our previous paper [TO] . For simplicity we will 
omit q A in the expression for U in (|1.3p . Taking into consideration the equalities 

S[Va( s a) = <(4) = K(sa\a, sa), S^Ksa) = s x ^° A (s A ), 
we obtain 

tyoAr^A = J>£Ka, sa) - e-^ w ^m{s k ))e-^ u ^ 

SA 
SA 

= ]T(e-^A) _ e -^(^))<( SA ) = 0. 

SA 

Here we changed signs of the spin variables sa in the first term in the sum in sa- 
Positive definiteness of the Hamiltonian follows from the following proposition. 

Proposition 3. The operator —Pa is positive definite on D(V A ). 

Proof. V A is an operator of multiplication by infinite differentiable functions on each one- 
dimensional spin subspace and its domain contains (<8>C 2 )' A ' ® Cq°(]R' a '). This domain coincides 
with the direct sum of 2l A l copies of of L 2 (lRl A l, e^°dq A ), where \q\o = \q x \\4> x \- 

The scalar product in the Hilbert space L A is given by (Fi,^) = J(Fi(q A ), F2(q A )) dq A , 
where the integration is performed over IRl A l We have to show that — (PAF(q A ), F(qa))o > 0- 
Let us define the operator 

p + = e lu{sl)p Ae -\u(sl)_ 
It is not difficult to check on the basis tp^ that 

where / is the unit operator. Here we used the following equality 

For the operator Pf we have 

P+F = J2(pXf)(<ia;sa)K(sa) 

SA 

and 

(P+F)(q A ;s A ) = -e-^M(F(q A ;s A )-F(q A ;si)). 
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It is convenient to introduce the new scalar product 

(Fi,F 2 )u = (e- u ^F 1 ,F 2 ) = X) / MQA,SA)F 2 (q A ,s A )e- U ^dqA 

SA J 

= J {F 1 {q A ),F 2 {q A )f u dq A = j \e~ u ^ F l (q A ), F 2 (q A )) dq A . 

The operator P^ is symmetric with respect to the new scalar product since 

(P+F 1 (q A ),F 2 (q A ))° u = (P A e-^ S ^F 1 (q A ),e-h u ( s X)F 2 (q A )) . (3.1) 

It is not difficult to check that 

-(P+F(q A ),F(qM = £ e^ u ^+ u ^ (F(q A ; s A ) - F(q A ; s£))F(q A ; s A ) 

SA 

= ^e-tM'Wrt^FfasJ - F(q A ;s£)) 2 > 0. (3.2) 

SA 

Here we took into account that the function under the sign of exponent is invariant under 
changing signs of the spin variables sa- Prom (|3.ip . (|3.2p it follows that —Pa > 0. I 

From Proposition [3] it follows that H A = JaPa is positive definite since for the nega- 

ACA 

tive J a we have 

(H{F,F) = J2 [ JA(qA)(PAF(q A ),F(q A )) dq A > 0. 
AC A 

The fact that the Hamiltonian is essentially self-adjoint is derived from the following propo- 
sition (see example X.9.3 in |18j). 

Proposition 4. The operator —A + fi 2 x 2 + V + (y,x), where (y,x) is the Euclidean scalar 

n 

product of x = Xj, j = 1, . . . , n with the constant vector y, A = Yl g§~> A* 7^ i/ y 7^ 0, V is 

i=i J 

the operator of multiplication by a non-negative function V(x) £ L 2 {W l ,e~ x dx), x 2 = (x,x), 
is essentially self-adjoint on S(M n ) PI D(V), i.e. the intersection of the Schwartz space and the 
domain of V . 

Proof. The following two inequalities are valid 

/i 2 ||x 2 Vll < ||(-A + y + / uV)V'll + 2/, 2 n||V||, (y,x)<ex 2 + ^-\y\ 2 , (3.3) 

4e 

where \i is a real number, \y\ = max|yL, || • || is the scalar product in the Hilbert space of 

j 

square integrable functions. We tacitly assume that x 2 means the operator of multiplication of 
a squared variable. The first inequality follows from the inequalities (dj = ^-) 

(_ A + V + fi 2 x 2 ) 2 > GuV) 2 - M 2 (Ax 2 + x 2 A), 
-(d)x\ + x\d 2 ) = -(2x k d 2 x k + 25 jM ) > -2S j>k , 

n 

-(Ax 2 + x 2 A) > - J^s? + x ffl)- 
3=1 



8 



T.C. Dorlas and W.I. Skrypnik 



Here we took into account that V is positive and [dj,Xk] = djXk — x^dj = 5j t k- As a result 

12 ' 

47 



\\(y,x)^\\<a\\(-A + V + ^x 2 )^\\+bU\\, a = ^" 2 e, 6 = n(2+ l2 



For \i 7^ number a can be arbitrary small and from the Kato-Rellich theorem [18] it follows 
that the essential domain of —A + V + fi 2 x 2 + (y, x) coincides with the essential domain of 
-A+V+fi 2 x 2 . But the last one coincides with §(R n )nD(l/) [El Theorem X.59]. Inequality (|3H 
is sufficient, also, for the proof of the proposition in the case fi = y = (this is explained in 
Example X.3 in [E])- ■ 

Since S A is a diagonal operator on ((X>C 2 )' A ' the operator 
- £ ^ + ^ 9 2 + V K + 2 W 2 £ &^(5i) 

sGA zGA xGA 

is the direct sum of 2l A l copies of the minus |A|-dimensional Laplacian plus the three functions 
coinciding with the three functions in Proposition [H From Proposition [J] it follows that this 
operator is essentially self-adjoint on the set ((g>C 2 ) |A| ®§(lRl A l) C[D(V\) that contains (<gC 2 ) |A| ® 
C§°(M' a '). The same is true for the Hamiltonian since operator (r///) 2 ^ fixi^A.) and the operator 

xeA 

depending on S\ in its expression are bounded. 

The operators h A , H\ are positive definite on the dense set ((X>C 2 )' A <8>S(IR' A ') f]D(V\) which 
is the essential set for H\. This implies that H A is positive definite on its domain D{H\) and >5 a 
is its ground state. 

Proof of uniqueness. We have to establish that the symmetric semigroup Pi, generated 
by —H\, maps non-negative functions into (strictly) positive functions (increases positivity) and 
this will imply that the ground state is unique (see Theorem XIII. 44 in [E])- We will establish 
this property with the help of a perturbation expansion. The kernel of the semigroup P\, 
generated by h\ + V\, is expressed in terms of the Feynman-Kac (FK) formula [EJCEH] 

where V^(q A ; s A ) = Va^a; sa) + Yl [m 2 0?x+#x(sa)) 2 -/•*]> P L q > (dw\) = l\ P\ xql (dw x ) is the 

xGA ' A x£A ' x 

conditional Wiener measure and w A (t) is the sequence of continuous paths. The semigroup P t 
is represented as a perturbation series in powers of Vq 

V o = ~Y1 JaS Uv 

ACA 

This series is convergent in the uniform operator norm [20] since Vq is a bounded operator. Its 
perturbation expansion is given by 



„ n+l 

P l = Y,Pl Pn= J dT 1 ---dr n P?l[(V P?- T *- 1 ), 



where r n +i = t . We now use the following simple inequality 



\Vj^(w\(T),s\)\dr 
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< V(wa) = \J\a 









r 


exp < 


aU + 









exp < aU + 2r]fi^2 4>x\w x {t)\ > + ^2(^ 2 (\w x (t)\ + r]<f) x ) 2 - //) 



xgA 



where U = max£/o(s A ), \J\\ = sup Y \ j a\ , <fe = max \4> x (sa)\. Let 
SA 9A aca SA 



xeA 



then it follows from the positivity of the the kernel P\ and Vo that 
P\g\, sa; q' A , «a) > A'(9a; g A ) V ^V?(* A ; 4)- 

Z * « .1 



n>0 



(3.4) 



Here we utilized the semigroup property of P t , the inequalities e a > e a ', Vo > V_, 

(VoP^-'KqA, s A ; q'a, s' A ) > P^- 1 (g A ; q'A)V-(s A ; , s' A ). 

Now, it is easily proved as in [10] that the matrix V_ is irreducible. As a result there exists 
a positive integer n such that and that V™, has positive non-diagonal elements [211 122j . Hence 
the kernel in the left-hand side of (|3.4p is positive. 



4 Order parameters 

In this section we will prove Theorem[Tl i.e. occurrence of different types of ho in the considered 
quantum systems. In the proof we will rely on the following basic theorem. 

Theorem 2. Let the ferromagnetic short-range potential energy U of the classical Ising model 
on the hyper-cubic lattice Z rf with the partition function Za = Y ex P{ — PU(sa)} be given by 

U(sA) = -Y,<P( A ) a lA]> 

ACA 

where <p(A) >0,s x = ±1. Let, also, the uniform bound (p(x,y) > (p > for nearest neighbors 
x, y hold and <p(A) = for A with odd number of sites. Then for a sufficiently large (p(3 > 1 
and the dimension d > 2 there is the ferromagnetic Iro, that is, for the Gibbsian two point spin 
average the uniform in A bound holds 

{o- x a y )A > 0, 

where ct x (sa) = s x and the magnetization (an order parameter) Ma = |A| _1 Y s x is non-zero 

xeA 

in the thermodynamic limit A — > Z d . 

The proof of this theorem is based on an application of the generalized Peierls principle 
(argument). It will be given in the end of this section (see, also, [171 126j . The next theorem is 
the consequence of the basic theorem. 

Proof of item I of Theorem Q3 Condition (|1.2|) shows that 

5^(a A ) = |A|+2£ MA)s [A] + J2\ Jo(x;A)s [A] 

x&A ACA xeA VzgACA 
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where J 2 (x U A) = Jo(x; A) and J 2 {A) = for odd \A\. The last term is equal to 
Yl J o( x ' A ) + 2 Yl Yl Jo(x;A 1 )J (x;A 2 )s [AlAA2] , 

xeAxgACA xeAx<£Ai,A2QA 

where = 1 and A\AA 2 = (A\ U A 2 )\(A\ D A 2 ). Due to translation invariance of interaction 

the first term is bounded by |A| ^ Jq(0;A), where the summation is performed over Z d , and 

A 

this expression is finite since the interaction is short-range. Hence subtracting a finite constant 
proportional to |A| from [/* one sees that the result admits representation (|4.1|) with positive J* 
instead of (p such that J* (A) = for odd \A\ and for nearest neighbors x, y the inequality 
J*{ x >y) > {2rj 2 n + a) J. The basic theorem and Proposition [1] imply occurrence of ferromagnetic 
ho in S* 3 and oscillator ho. ■ 

Proof of item II of Theorem [TJ, Let iPa{qa] s\) = ^a( s a)V , oa(?a)- Then (|1.2|) and the defi- 
nition of S" 1 imply that 

»A «A 

Taking into account also the orthonormality of the basis one obtains 

2 

= ^-!y^ e t '6A e -f [^o(s A ' y )+^o(sA)] 

2 , 



*A 



e -aU {s A ) _ 



«A 

We also have 

^(^(^A) + ^(4' y )) 2 

x'gA 

= £ [3^„( SA ) + ^,( a JV)] + 2 £ M*a)["M«a) + 
x'eA x'eA 



a-'GA x'GA x'GA 



>4^^,( Sa )-B 2 -2C7Bi. 



of/ ( SA ) 



x'GA 

This yields 

{SlS l y ) A > e-W^e-^ZX 1 £ e * V 

«A 

= e -4 ii (- B 2+ 2CB i)e-f Bo . I 
Proof of Proposition [21 It is obvious that |^>x(s A )| < ||^o||i an d that 

- <M«a)| = I - 2[s y J (y - x') + J (x - x')]\ < 2[\J (y - x')\ + \J (x - x')\]. 
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As a result Wx]y(s\) < 4||J ||i. Further 

|^( S ^)-^ ; ( SA )| = 



>] M z - x ') s z - s y J o(y - x ') - s x Jo(x - x') 

zeA\(x,y) 

- 2 

Jq{z - x')s z + s y J (y - x') + s x J (x - x') 

zeA\(x,y) 

-4 ^2 M z ~ x ') s z(s y J {y - x') + s x J (x - x')) 

zeA\(x,y) 

< 4^|J (z - x')\(\J (y -x')\ + | J (a - x')\) < 4|| J ||i(| J (y -x')\ + \J (x - x')\). 

zeA 

Hence wi2(s A ) <8| I Jo 1 1?. ■ 
Proof of Theorem [2j Let xt = \(\ ^ a x) then one obtains 

MxtXy )A = 1 + (^)A - K)A - (o- x (Ty)A- 

Since the systems are invariant under the transformation of changing signs of spins the third 
and the second terms in the right-hand side of last equality are equal to zero and 

{<r x <ry)A = 1-4(x^X«)a- 



Hence if 

(xtXy }A < ^ 

then the ferromagnetic ho occurs, i.e. 
(Ox<?y)A > a > 0, 



(4.2) 



(4.3) 



where a is independent of A. If one succeeds in proving that there exists a positive function Eq((3) 
and positive constants a, a' independent of A such that 



(xtXyh < a'e E ^ 



(4.4) 



and proves that Eq is increasing at infinity then (|4.2p . (|4.3p will hold for a sufficiently large 
inverse temperature f3. The Peierls principle reduces the derivation of (|4.4p to the derivation of 
the contour bound. ■ 



Peierls principle. Let the contour bound hold 



(4.5) 



<x )3/ >er 



where (-)a denotes the Gibbs average for the spin system confined to a compact domain A, 
r is a set of the nearest neighbors, adjacent to the (connected) contour, i.e. a boundary of the 
connected set of unit hypercubes centered at lattice sites. Then (|4.4|) is valid with Eq = a"E, 
where a" is a positive constant independent of A. 
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Proof of contour bound. Bricmont and Fontain derived the contour bound for the spin sys- 
tems with the potential energy (|4.1|) with the help of the second Griffiths [23] and Jensen 
inequalities [23] (see also [251 12S] ) 

(^[A]^[b])a[t] - (o>])A[r](o-[i?]>A[r] > 0, j e f dfi > exp j ^ fd/A , 

where d[i is a probability measure on a measurable space. Their proof starts form the inequality 

xtXy = e-i^ef^x+x- < e'^'^xtXy < <r*°* ay ■ 
As a result (/?' = (p[3) 

( II xtXyh < (e <^ >er )a = (e <-> er )X[ r] 
(x,y)er 

~4~ S ( cr a: cr y>A[r] 

<e {x - y}er =e~ Ev , 



where (•, -)A[r] is the average corresponding to the potential energy 
Uv{qh) = U(s A ) + — ^2 SxS v 

In the last line we applied the Jensen inequality. From the second Griffiths inequality it follows 
that the average (cr x & y )A[r] is a monotone increasing function in ipA- So, in the potential energy 
determining this average we can put ipA = 0, except A = (x, y) and leave the coefficient (p in 
front of the bilinear nearest-neighbor pair potential in ()4.ip without increasing the average. This 
leads to 

{a x a y ) m >{aa') = Z^(^\ ]T SlS2 e^ Z 2 ((3) = £ e 4r*i-». 

S1,S2=±1 Si,S2=±l 

That is, 

E r > \T\E, E = 2~ 1 p'{aa') 



or 

E = (5{e 2 -^' - e- 2 " 1 0(e 2 " 1 ' 9 ' + e" 2 " 1 ' 3 ')- 1 > 2~ 1 /3'(1 - e^'). 

Here we used in the denominator the inequality e~ 2 < e 2 & . Obviously, E tends to infinity 
if 0' tends to infinity. This implies (|4.5j) . ■ 



5 Discussion 

We showed that in the considered lattice spin-boson models with J a < ground states are 
Gibbsian and the ground state averages for special observables are reduced to averages in classical 
Ising models. This means that existence of ground states order parameters is connected with 
existence of order parameters in the associated Ising models and that a breakdown of symmetries 
in the quantum systems is determined by a breakdown of symmetries in Ising models. We 
considered the free boundary conditions implying that for the cases of the perturbation V\, 
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considered in the two theorems, the ground state averages of q x , S% are zero, that is (q x )k = 0, 
(S^)a = if the associated Ising potential energy is an even function. In order to make such the 
averages non-zero (explicit symmetry breaking) one has to introduce special boundary conditions 
(quasi- averages) which have to single out pure Gibbsian states in the associated Ising models. It 
is known [28] that for the two-dimensional ferromagnetic Ising nearest-neighbor model there are 
two boundary conditions which generate pure states and that every other state is a convex linear 
combination of these two states. A discussion of a construction of ground states in lattice spin 
and fermion quantum systems with an explicit symmetry breaking a reader may find in |29j . 

Remark 2. Translation invariance means that 

Jxi,...,x n — Jo,X2—xi,...,x n —xi j JQ\pC'i %1 j ■ ■ ■ > %n) — ^o(O) ^1 X, . . . , X n x) 

where J, Jo are symmetric functions. The short-range character of interaction means that 

max \Jx,a\ < QO; max |Jp(x; A)\ < oo. 
A A 

Remark 3. If only one-point sets are left in the sum for V\ then the expression for H\ can be 
rewritten in the following way 



The property of the ground state to be a ground state with the zero eigenvalue of a local 
Hamiltonian H x was found earlier for special isotropic anti-ferromagnetic Heisenberg chains with 
valence bond ground state in |27j . 
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